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1 Introduction 



Lax representation of classical dynamical evolution equations [T] is one key ingredient in 
the modern theory of classical integrable systems [2]^[Z1 together with the associated notion 
of classical r- matrix [H [9] . It takes the generic form of an isospectral evolution equation: 
^ = [L, A], where L encapsulates the dynamical variables and A defines the time evo- 
lution. We shall generically consider the situation where L and A depend on a complex 
(spectral) parameter. The spectrum of the Lax matrix or its extension (transfer matrices), 
or equivalently the invariant coefficients of the characteristic determinant, thus provide au- 
tomatically candidates to realize the hierarchy of Poisson-commuting Hamiltonians required 
by Liouville's theorem [TT]. Existence of the classical r-matrix then guarantees Poisson- 
commutativity of these natural dynamical quantities taken as generators of the algebra of 
classical conserved charges. 

The question of finding an appropriate time evolution matrix A for a given Lax matrix L 
therefore entails the possibility of systematically constructing classically integrable models 
once a Lax matrix is defined with suitable properties. Depending on which properties are 
emphasized, this problem may be approached in several ways. One approach uses postulated 
Poisson algebra properties of L, specifically the r-matrix structure, as a starting point, and 
establishes a systematic construction of the time evolution operator A associated to the 
Hamiltonian evolution obtained from any function in the enveloping algebra of the Poisson- 
commuting traces of such an L matrix. Such a formulation was proposed long time ago [HI [9] 
for the bulk case, when the Poisson structure for L is a simple linear or quadratic r-matrix 
structure. 

We shall consider here more general situations when the dynamical system also depends 
on supplementary parameters encapsulated into a matrix K. We shall restrict ourselves 
to the situation where these parameters are non-dynamical, i.e. the Poisson brackets with 
themselves and with the initial "bulk" parameters is zero. For situations where the extra pa- 
rameters are dynamical see [T2l[T3] . Note that interpretation of the physical meaning of these 
c-number parameters will come a posteriori when computing the associated Hamiltonians. 
In particular any physical interpretation of the K matrix as a description of the "boundary 
properties" (external fields, ...) may not be appropriate in all cases as shall appear in our 
discussion of examples. We shall however keep this designation as a book-keeping device 
throughout this paper. 

Our central purpose will be twofold. We shall first of all define generic sets of sufficient 
algebraic conditions on K also formulated in terms of the "bulk" r-matrix. A "boundary- 
modified" generating matrix of candidate conserved quantities, hereafter denoted T, will be 
accordingly constructed as a suitable combination of L and K matrices. The idea for such 
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a construction naturally arises when considering the semi-classical limit of the well-known 
Cherednik-Sklyanin reflection algebras preserving bulk quantum integrability [11]. We shall 
then redefine the time evolution operator A associated to a given Hamiltonian constructed 
from T (modified monodromy matrix) from the new basic elements, i.e. the matrix T, 
the bulk r- matrix or (r, s) pair [9l [15], and the reflection matrix K. Such a construction 
was exemplified in [TB]; what we propose here is however a generic, systematic procedure 
to obtain modified "boundary" or "folded" classical integrable systems from initial pure 
"bulk" systems. Note also that the example worked out in [16] is precisely a case which the 
X-matrix does describe boundary effects. 

It must be emphasized here that an alternative, analytic approach to this question, at 
least in the bulk case, was extensively described in [7] (chapter 3). In this approach one 
uses instead the analytic properties (location of poles and algebraic structure of residues) 
of the Lax matrix L as a meromorphic function of the spectral parameter. They provide 
for a unique consistent form of an associated A matrix. The subsequent Lax equation is 
then developed into separate equations corresponding to the poles of L and A. The Poisson 
structure and r matrix structure are only then defined as providing a consistent Hamiltonian 
interpretation of this Lax equation. Similar reformulation of our results must exist but we 
shall not consider them here. 

We shall expand here the r-matrix approach in three situations. We start with the 
simpler case where the initial bulk structure is a linear Poisson structure for a Lax ma- 
trix parametrized by a single r-matrix. We then develop the cases of both discrete and 
continuous parametrized Lax matrices relevant to the description of systems on a lattice 
or on a continuous line. In these last two situations the relevant r-matrix structure is a 
quadratic Sklyanin-type bracket [17] . We shall in all three cases derive (or actually rederive, 
in some cases) the form of the generating functional for Poisson-commuting Hamiltonians, 
and establish the explicit general formula yielding the A operator. Explicit examples shall be 
developed in all cases, albeit restricting ourselves to the simplest situations of non-dynamical 
non-constant double-pole rational r-matrices. More complicated situations (trigonometric 
and/or dynamical r-matrices, relevant for e.g. sine-Gordon models or affine Toda field the- 
ories) will be left for further studies. 

2 Linear Poisson structure 

We consider here the original situation [1] where the full dynamical system under study is 
represented by a single Lax matrix, living in a representation of a finite-dimensional Lie alge- 
bra or a loop algebra. In this last case the Lax matrix also depends on a complex parameter 
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A known as "spectral parameter" . In all cases the requirement of Poisson-commutativity for 
candidate Hamiltonians TrL^ is equivalent [18] to the existence of a classical r-matrix [HI [9] 
realizing a linear Poisson structure for L. Indeed, consider the Lax pair (L, ^4) satisfying 

dL 
'dt 

the associated spectral problem 

L(A) if) = u ip, 



A, L 



2.1 



det(L(A) -u)=0 



(2.2) 



provides the integrals of motion, obtained through the expansion in powers of the spectral 
parameter A of trL{X). Alternatively, in particular when no spectral parameter exists, one 
should consider the traces of powers of the Lax matrix trL"(A) as natural Hamiltonians. 

Assuming that the L matrix satisfies the fundamental relation 



rabi^ - At), La{X) + Lbifl) 



it is shown using (12. 3p that for any integer n,m: 

|trL"(A), trL'"(^)} = 0. 



(2.3) 



(2.4) 



The reciprocal property was shown in [18] . Let us recall [H [9], [18] how one may identify the 
74-operator associated to the various charges in involution. Bearing in mind the fundamental 
relation (12.31) it is shown that 



tra L^(A), Lfe(^)} = n tva [l^-^X) rat{X - /i)) L,{^i) - n L,{fi) tva (^L:-\X)rat{X - /i))(2.5) 
From (12. 5p one extracts A: 

Ar,{X, fi) = ntra (^L:-\X)r,,{X-fi)y 
In the case of the simplest rational non-dynamical r-matrices [I9] 



(2.6) 



r(A) 



P 
A 



N 



where P = Eij 



E. 



(2.7) 



P is the permutation operator, and {Eij)ki = SikSji, we end up with a simple form for A^. 



n 



X — ji 



L^~\X) 



(2i 



and as usual to obtain the Lax pair associated to each local integral of motion one has to 



expand An = Note that generically, using the dual formulation of the classical 

r-matrix |9l one also has: 



A{X,^) = Tr{r{X, fi)dH) 



(2.9) 



3 



where H is the Hamiltonian expressed as any function in the enveloping algebra generated 
by TrL*^. 

Ultimately we would like to consider an extended classical algebra in analogy to the quan- 
tum boundary algebras arising in integrable systems with non-trivial boundary conditions 
that preserve integrability. Subsequently we shall deal with two types of algebras which may 
be associated with the two types of known quantum boundary conditions. These boundary 
conditions are known as soliton preserving (SP), traditionally studied in the framework of 
integrable quantum spin chains (see e.g. [H], [2D]-[22]), and soliton non-preserving (SNP) 
originally introduced in the context of classical integrable field theories [16], and further 
investigated in [2U[25]. SNP boundary conditions have been also introduced and studied for 
integrable quantum lattice systems [26]-[30]. From the algebraic perspective the two types 
of boundary conditions are associated with two distinct algebras, i.e. the reflection algebra 
[H] and the twisted Yangian respectively [3ll[32] (see also [25l[29l[30l[33l[3ll). The classical 
versions of both algebras will be defined subsequently in the text (see section 3.2). It will 
be convenient for our purposes here to introduce some useful notation: 

r,,(A) forSP, f,,(A) = r*f "(A) for SNP 
TabiX) forSP, r:,(A) = r*^„(-A) for SNP 
na{X) forSP, f:,(A) = r*",(-A) for SNP (2.10) 

together with: 

L(A) = -L(-A) for SP, L(A) = L*(-A) for SNP. (2.11) 

In addition the "boundary conditions" , to be interpreted on specific examples, are parametrized 
by a single non-dynamical matrix k(A). We propose here the following set of algebraic rela- 
tions: 

{ri(Ai), %{X2)}=T^2iXuX2) Ti(Ai)-Ti(Ai) fr2(Ai,A2) 

+r+ (Ai, As) r2(A2) - r2(A2) i^X,, X^) (2.12) 

where we define: 

ri2(Ai - A2) k2(A2) - k2(A2) r*2(Ai + A2) 
k2(A2) ri2(Ai - A2) - r*2(Ai + A2) k2(A2) 
ruiXi - A2) ki(Ai) + ki(Ai) f*2(Ai + A2) 
ki(Ai) fi2(Ai - A2) + rl^iXi + A2) ki(Ai), (2.13) 

for a k matrix satisfying: 

|ki(A), k^ifi)} = 0, {ki(A), L2(/i)} = 0, (2.14) 



rab{X) = 
r:,{X) = 



Ti2[Xi, A2) = 
fl2('^l' ^2) = 

f^2('^li ^^2) = 
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ri2(Ai - A2) ki(A) k2(A2) + ki(Ai) f*2(Ai + A2)k2(A2) 

= ki(Ai) k2(A2) fi2(Ai - A2) + k2(A2) rt2(Ai + A2) ki(Ai). (2.15) 

We then show: 

Theorem 2.I. : The quantity 

r(A) = L(A) k(A) + k(A) L(A) (2.16) 
is a representation of the algebra defined by (12.121) . with k obeying (12.141) . (I2.15p . 

Proof: We shall need in addition to (12. 3p the following exchange relations: 

L„(A), Lb(/i)} = L„(A)f:,(A + + f:,(A + - fl^X + /i)L,(A) - Lb{fi)r:,{X + /x) 

L,(A), L,(/i)} = L„(A)r:,(A + ^) + r:,(A + - r:,(A + /x)L,(A) - L,(/i)r:,(A + /i) 

'l„(A), L,(/i)| = kfc(A - /i), L,(A) + Lb(^)l . (2.17) 



By explicit use of the algebraic relations (12.151) and (12.171) we obtain: 

L„(A)k,(A) + L„(A)k,(A), U{fi)h{fi) + h{fi)Lb{fi) 
= r-,(A,/i)(L,(A)k„(A) +k,(A)L,(A)) - (L,(A)k,(A) + k„(A)L,(A) jf", (A, /i) 
+ra6(A,/i)(^b(^)k6(^) +k6(/i)L6(^)j - (^Lb(^)k6(^) +k6(/i)L6(^)jf+(A,/i) 
+ (^LaiX) + Lbifi)^ ( - rabiX - iJ,)kaiX)kb{fi) - ka(A)f*b(A + /i)kb(/i) 
+kb{fi)rl^{X + /i)k„(A) + ka(A)kfe(/i)fi2(A - /i)) 
+ ( - rab{X - fi)ka{X)kb{fi) - ka(A)f*^(A + n)kb{fi) 

+kb{fi)rl^{X + /x)k,(A) + k,(A)k,(/x)fi2(A - /i)) (l„(A) + L,(/i)) . (2.18) 

Bearing however in mind that the k-matrix obeys (12.151) we conclude that the last four lines 
of the equation above disappear, which shows that (I2.16P satisfies (I2.12p . and this concludes 
our proof. ■ 

Define T(A) = k~^(A) T(A), we now prove: 
Theorem 2.2.: 

'trX(A), tr,Tf(/.)|=0. (2.19) 
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Proof: 

employing fl2.12p the preceding expression becomes 
. . . oc 

= ... = 0. (2.21 

Note that in order to show that the latter expression is zero we moved appropriately the 
factors in the products within the trace and we used (12. 151] . ■ 

We finally identify the modified Lax formulation associated to the generalized algebra 
(^Ml-I^J^ as: 

Theorem 2.3. : Defining Hamiltonians as: tr(jT"(A) = the classical equations of mo- 

tion for T: 

= T(/x)} (2.22) 

take a zero curvature form 

t(/i)= [a(A,aO, T(/i)], (2.23) 

where A„ is identified as: 

K{\lA = n tr,(Tri(A)k:i(A)f+(A,^)). (2.24) 

Proof. 

'ir,T:(A), T,(^) 



n tr„(Tri(A)k;i(A)k-i(/i)(r-(A,^)r.(A) -r„(A)f-(A,^) + r+(A,^)r,(/i) -r,(/i)f+ 

= n tr,(Trik-i(A)k,(/i)(r;,(A,^)T„(A) - ,,))) 

+n tr„(Trik:i(A)k,(^)(r+(A,/.)T,(^) -T,(/i)f+(A,/x))). (2.25) 

Taking into account (12.151) we see that the first term of RHS of the equality above disappears 
and the last term may be appropriately rewritten such as: 

tr,T:(A), T,(^)} = n tr, (T"-i(A)k-i(A)f+ (A, ^)) T,(/i) - n T,(/i) tr, (T"-i(A)k,-i(A)f+ (A, ^)) . 



(2.26) 
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From the latter formula fl2.24p is deduced. ■ 

Finally, let T(A), T'(A) be two representations of (12.121) . and let also 

{r„(A), TM}=0- (2.27) 

It is then straightforward to show, based solely on the fact that T, T' satisfy (I2.12p and 
(I2.27P , that the sum T(A) + T'(A) is also a representation of (I2.12p . 



2.1 Examples 

We shall present here a simple example, starting from the classical rational Gaudin model 
|35j . Details on the so called "dual" description of the Toda chain [361 EZ] and the DST 
model [HHl [39], both associated to the A^^_^ r-matrix |40], will be presented in a forth- 
coming publication. Before we proceed with the particular example let us rewrite A„. We 
consider a situation where the Poisson brackets are parametrized by the simplest rational 
non-dynamical r-matrix [19]. After substituting the rational r-matrix in (12.241) we get for 
both types of boundary conditions: 

V A — yU A + /i 

The L-matrix associated to the classical glj\/ Gaudin model, and satisfying the algebraic 
relation (12.31) is 

^(^) = E E X^.^-^ (2.29) 

a,/3=l 1=1 

where Vai3 G gW- Recall that in the "bulk" case the integrals of motion are obtained through 
tr L^i^X). The first non-trivial example reads 

t(^)(A) = tr L\X) =J:J: (2.30) 

and by taking the residue of the latter expression for A — we obtain the Gaudin Hamil- 
tonian: 

=Y.Y. -fzf- (2.31) 
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Let us now come to the generic algebra (12.121) . (12.141) . (I2.15P considering both SP and SNP 
cases. Based on the definitions (12.111) we have: 

M N Q(i) 

= EEaT7^« -here 

a,l3=l i=l * 

Q.(5 = V^p forSP, Qo.(, = -Vp^ for SNP. (2.32) 

Let us take as a representation T(A) = L(A)k(A) + k(A)L(A) + K{\) where k is a solution 
of (I2.14p . (I2.15p . and i^' is a c-number representation of (I2.12p with zero Poisson bracket. 
To obtain the relevant Hamiltonian we now formulate traT^(A) (where T = k~^T and we 
also define K = k~^K), and the Hamiltonian arises as the residue of the latter expression at 
A = zy. 



i^j=la,f3=l * ^ j ,i=l a,/3,7,5,e=l * ^ i=l a,l3,y,S,e=l 



A special case of the generic algebra (12.121) is discussed in [ITl W2\ . Note also that expression 
(I2.33P may be also obtained as the semiclassical limit of the quantum gl^ inhomogeneous 
open spin chain for special boundary conditions, (see e.g. [43]), zj being the inhomogeneities 
attached to each site j. However, here the Hamiltonians are obtained directly from the 
representations of our new classical algebra (12.121) . (12.141) . (I2.15p . It appears in this example 
that the parameters k(A) play the role of "coupling constants" consistent with the integrable 
"folding" of a 2N site Gaudin model, and not the role of boundary parameters. 



3 Quadratic Poisson structures: the discrete case 

Quadratic Poisson structures first appeared as the well-known Sklyanin bracket [Hj. A 
more general form, characterized by a pair of respectively skew symmetric and symmetric 
matrices (r, s) appeared in [15] in the formulation of consistent Poisson structures for non- 
ultralocal classical integrable field theories. Finally it was shown ^5] that this was the 
natural quadratic form a la Sklyanin for a non-skew-symmetric r-matrix, reading: 

L2} = [r - r", L1L2] + Li(r + r'^)L2 - L^ir + r^)L,. (3.1) 

A typical situation when one considers naturally a quadratic Poisson structure for the Lax 
matrix occurs when considering discrete (on a lattice) or continuous (on a line) integrable 
systems where the Lax matrix depends on either a discrete or a continuous variable; the Lax 
pair is thus associated to a point on the space- like lattice or continuous line [6l [IT]. Let us 
first examine the discrete case where one considers a finite set of Lax matrices L„ labelled 
by n G N. 
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3.1 Periodic boundary conditions 



Introduce the Lax pair (L, ^4) for discrete integrable models [1] (see also [IB] for statistical 
systems), and the associated auxiliary problem (see e.g. [6]) 



(3.2) 



From the latter equations one may immediately obtain the discrete zero curvature condition: 



The monodromy matrix arises from the first equation (13.21) (see e.g. [6]) 

T,(A) =L,;v(A)...L,i(A) 



(3.3) 



(3.4) 



where index a denotes the auxiliary space, and the indices 1, . . . , denote the sites of the 
one dimensional classical discrete model. 

Consider now a skew symmetric classical r-matrix which is a solution of the classical 
Yang-Baxter equation P E] 



ri2(Ai - As), ri3(Ai) + r23(A2) + ri3(Ai), r23(A2) 
and let L satisfy the associated Sklyanin bracket 

La(A), Lfc(^)[= rab(A-/i), La{X)Lb{iJ,) 



0, (3.5) 



(3.6) 

It is then immediate that that (13.41) also satisfies (13. 6p . Use of the latter equation shows that 
the quantities trT(A)" provide charges in involution, that is 



tr T"(A), tr T'"(/x)} =0 



(3.7) 



which again is trivial by virtue of (13.61) . In the simple s/2 case the only non trivial quantity 
is trT{X) = t{X), that is the usual "bulk" transfer matrix. 

Let us now briefly review how the classical A-operator and the corresponding classical 
equations of motion are obtained from the expansion of the monodromy matrix in the simple 
case of periodic boundary conditions. Let us first introduce some useful notation 



Ta{n, m; A) = La„(A)La„_i(A) . . . ^^^(A), n > m. 



{31 



To formulate {t(A), L{fi)} or even better {In t{X), L(/i)}, given that usually In trT{X) gives 
rise to local integrals of motion, we first derive the quantity {Ta{X), Lbn{fi)}. 
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|Ta(A), Lf,„(//)| = Ta{N,n + l;X) TabiX- n) Ta{n,l;X) LbnilJ.) 

- LM Ta{N,n;X)rabiX-fi) Tain -1,1;X). (3.9) 



Taking the trace over the auxihary space a, and then the logarithm we conclude 

{in t{X), L(^)} = t-\X) (tra{Ta{N,n + l;X) rab{X-^i) T,(n,l;A)} L(//) 
- L{iJi) tr,{Ta{N,n- A) r,b{X - ^i) Ta{n - 1, 1; A)}' 



(3.10) 

the auxiliary index h is suppressed from the latter expression for simplicity. Then plugging 
the latter forms into the classical equations of motions for all integrals of motion 

LM = {t(A), L„(^)} (3.11) 

we conclude that the zero curvature condition (13. 3p is realized by: 

An{X, ii) = r\X) tra{Ta{N, n- A) rab{X - ^l) Ta{n - 1, 1; A)}. (3.12) 

Let us focus on the classical rational r- matrix (12.71) . In this case An takes the simple form: 

An{X, fi) = Tin - 1, 1; A) TiN, n; A). (3.13) 
A — /i 

In the open case, as we shall see in the subsequent section, it is sufficient to consider the 
expansion of t{X) in order to obtain the local integrals of motion. 



3.2 Open boundary conditions 

We now generalize the procedure described in the preceding section to the case of generic 
integrable "boundary conditions". We propose a construction of two types of monodromy 
and transfer matrices, and associated Lax-type evolution equations, of the form (I3.6p - (l3.12p . 
albeit incorporating a supplementary set of non-dynamical parameters encapsulated into a 
"reflection" matrix K{X). In some examples they may indeed be interpreted as boundary 
effects consistent with integrability of an open chain-like system. We should stress that this 
is the flrst time to our knowledge that such an investigation is systematically undertaken. 
There are related studies regarding particular examples of open spin chains such as XXZ, 
XYZ and ID Hubbard models [171 HE]. However, the derivation of the corresponding Lax 
pair is restricted to the Hamiltonian only and not to all associated integrals of motion of 
the open chain. In this study we present a generic description independent of the choice of 
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model, and we derive the Lax pair for each one of the entailed boundary integrals of motion. 
Particular examples are also presented. 

The two types of monodromy matrices will respectively represent the classical version of 
the reflection algebra M, and the twisted Yangian T written in the compact form: (see e.g. 

{ri(Ai), T2(A2)} = ri2(Ai - A2)Ti(Ai)r2(A2) - ri(Ai)T2(A2)fi2(Ai - A2) 
+Ti(Ai)f*2(Ai + A2)T2(A2) - T2(A2)r*2(Ai + A2)Ti(Ai) (3.14) 

where f, r*, r* are defined in (12.101) . The latter equation may be thought of as the semi- 
classical limit of the reflection equation [M] . In most well known physical cases, such as the 
^^f-l r-matrices rjg*^ = r2i implying that in the SNP case r*^ = f*^. In the case of the 
Yangian r-matrix ri2 = r2i, hence all the expressions above may be written in a more sym- 
metric form. These two distinct algebras are respectively associated with the SP boundary 
conditions (M) and the SNP boundary conditions (T). 

In order to construct representations of (13.141) yielding the generating function of the 
integrals of motion one now introduces c-number representations (K^) of the algebra M (T) 
satisfying (13.141) for SP and SNP respectively, and also the non-dynamical condition: 

{KtiXi), K^{X2)}=0. (3.15) 

Taking now as T(A) any bulk monodromy matrix (13.41) built from local L matrices obeying 
(13. 6p and defining in addition 



f{X)=T~\-X) forSP, r(A)=T*(-A) for SNP. (3.16) 

one gets: 

Theorem 3.I. : Representations of the corresponding algebras R, T, are given by the fol- 
lowing expression see e.g. [T^ [13]: 



r(A) = r(A) K'{X) T(A). (3.17) 

For a detailed proof see e.g. [13]. ■ 

Define now as generating function of the involutive quantities 

t{X) =tr{K+{X)T{X)}. (3.18) 



Theorem 3.2. : Due to fIXTD it is shown that [H, [13] 

|t(Ai), t(A2)} = 0, Ai, A2GC. 
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(3.19) 



The expansion of t{X) naturally gives rise to the integrals of motions i.e. 

^(^) = E^- (3-20) 

i 

Usually one considers the quantity In t(A) to get local integrals of motion, however for the 
examples we are going to examine here the expansion of t(A) is enough to provide the asso- 
ciated local quantities as will be transparent in the subsequent section. Finally one has: 

Theorem 3.3. : Time evolution of the local Lax matrix L„ under generating Hamiltonian 
action of t{\) is given by: 

L„(/i) = A„+i(A, ^) Lniii) - A„(A, yu) (3.21) 

where A„ is the modified (boundary) quantity, 

A„(A, ^) = tra{Kt{\) Ta{N, n; A) r,b(A - /i) r,(n - 1, 1; A) /C(A) fa{\) 

+ Kt{\) U\) K-{\) t(l,n-l;A) f:,(A + /i) f;(n,iV;A)). (3.22) 

Proof. We need in addition to (13. 6p one more fundamental relation i.e. 

{4(A), L,(^)} = L,(A)f:,(A)L,(^) -L,(^)f:,(A + /x)4(A). (3.23) 

where the notation L is self-explanatory from (13.161) . Taking into account the latter expres- 
sions we derive for T 



Ta{X), LM\ = T„(l,n;A) f:,(A + /x) T,(ri + l,iV;A) L,„(/i) 

- Lfc„(/i) f,(l,n-l;A) f:,(A + /i) T;(n,iV;A). (3.24) 



The next step is to formulate |t(A), L6„(/i)|; indeed recalling ([33]), ([321, ff^TH) and fl3A8l) 
we conclude: 

\t{\), LM] = tr,(/C(A) T,(iV,n + l;A) r„b(A-/i) T,(n,l;A) i^-(A) r(A) 

+ /C(A) T,(A) JC(A)Ta(l,n;A) f:,(A + ^) t(n + 1, iV; A))Lfc„(/i) 
- LM tra[Kt{X) T,iN,n; A) r,,{X - fi) T,{n - 1, 1; A) i^,^(A)r(A) 
+ K+(A) T„(A) K-{\) t(l,n- 1;A) f:,(A + ^) T;(n,iV;A)). (3.25) 
Expression (13.221) is readily extracted from (13.251) . ■ 
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Special care should be taken at the boundary points n = 1 and n = N + 1. Indeed going 
back to formula (13.221) restricting ourselves to n = 1 and n = + 1 and taking into account 
that T(A^, A) = T(0, 1,A) = r(l,0,A) =r(A^ + l,A^, A) = I we obtain the explicit form 

for Ai, A„+i. We should stress that the derivation of the boundary Lax pair is universal, 
namely the expressions (13.221) are generic and independent of the choice of L, r. Note that 
expansion of the latter expressions of A„(A,//) = ^ ^ provides all the A^*'' associated to 
the corresponding integrals of motion 7i^*\ which follow from the expansion of t(A). This 
will become quite transparent in the examples presented in the subsequent section. 

Remark: A different construction of (13.141) was already given in a very general setting 
in [l9] . It is related to the formulation of non-ultralocal integrable field theories on a lattice 
and extends the analysis of [H]. The essential difference with our construction is that the k 
matrix (denoted 7) in [19], is sandwiched between the "local" monodromy matrices Tn.n-i 
so as to obtain an overall Poisson bracket of the form (I3.14p for the "dressed" monodromy 
matrix ...Tn+i^njTn^n-il---- The matrix 7 allows to take into account the effects of the non- 
ultralocal part 6'{x — y) of the Poisson bracket structure. In this respect 7 also corresponds 
to "boundary" effects, although multiple and internal. 

3.3 Examples 

We shall now examine a simple example, i.e. the open generalized DST model, which may be 
seen as a lattice version of the generalized (vector) NLS model, (see also [391 Eni EH [521 [13] 
for further details). The open Toda chain will be also discussed as a limit of the DST model. 
We shall explicitly evaluate the "boundary" Lax pairs for the first integrals of motion. If we 
focus on the special case (12.71) . which will be our main interest here, the latter expression 
reduces to the following expressions for SP and SNP boundary conditions. In particular, A„ 
for SP boundary conditions reads: 

AniX,fi) = T{n- 1,1; X) K-{X)f{X) K+{X)T{N,n;X) 

X — jj, 

+ fin, N; A) K^(X) T(X) K~(X) f(l, n - 1; A) for SP (3.26) 

A + /i 

and for SNP boundary conditions: 

A„(A,/i) = -^T{n-l,l;X) K"{X)f{X) K+{X)T{N,n;X) 

A — jJ, 

- -^f\l,n-l;X) K-\X)T\X) K^\X)f\n,N;X) for SNP, 
A + /i 

(3.27) 

where we recall that for the special points n = 1, 1 we take into account that: T{N, N + 
1, A) = r(0, 1, A) = f (1, 0, A) = f{N + 1, N, A) = I. 
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The Lax operator of the gl{M) DST model has the following form: 

AT-l MM M 

L(A) = (A - J] x(^)x(^))£;n + b^E^, + b^x^^~^'^E^j -^X^^-'^E^^ (3.28) 

j=l j=2 j=2 j=2 

with xii\ xii'^ being; canonical variables 

I, jG{l,...,Ar}, n, mG{l,...,iV}. (3.29) 

In [13] the first non-trivial integral of motion for the SNP case, choosing the simplest con- 
sistent value = I was explicitly computed: 

^ N N M-1 . M-1 

n=l n=l j=l i=l 

M~l 

where nn = J2^n^^n^- (3-30) 

Our aim is now to determine the modified Lax pair induced by the non-trivial integrable 
boundary conditions. We shall focus here on the case of SNP boundary conditions, basically 
because in the particular example we consider here such boundary conditions are technically 
easier to study. Moreover, the SNP boundary conditions have not been so much analyzed in 
the context of lattice integrable models, which provides an extra motivation to investigate 
them. Taking into account (13.271) we explicitly derive the modified Lax pair for the general- 
ized DST model with SNP boundary conditions. Indeed, after expanding (13.27P in powers 
of we obtain the quantity associated to the Hamiltonian (I3.30I] R: 

A(f) = Ai?n - E^-i'^^^i + bJ2^n''^Ey, ne{2,...N} 

^^N+i = ^-^11 - ^^-^Ji + E^^ ^^E^j- (3.31) 

Let us now consider the simplest possible case, i.e. the s/2 DST model. It is worth stressing 
that in this case the SP and SNP boundary coincide given that 

L~\-X) = V L\-X) V, \/ = antid(l,...,l). (3.32) 



■^To obtain both Ti*^^^ and aI^'' we divided the original expressions by a factor two. 
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In this particular case the Hamiltonian fl3.30p reduces into: 

1 ^ }? 1 

n=l n=l 

The equations of motion associated to the latter Hamiltonian may be readily extracted by 
virtue oi 



L = 

It is deduced from fl3.3ip for the s/2 case: 



n^^\ l}. (3.35) 



Alternatively the equations of motion may be derived from the zero curvature condition 

^ = A« , L„ - L„A» (3.37) 

which the modified Lax pair satisfies. It is clear that to each one of the higher local charges 
a different quantify An^ is associated. Both equations (13.351) . (13.371) lead naturally to the 
same equations of motion, which for this particular example read as: 

Xn = xlX„, + bXn+1, X„, = -XnXl - bX^-i, u G {2, . . . N - 1} 
Xi = x\Xi + hx2, Xi = —XiXl — hxi 

xtv = x\X]^ + X]s[^ X]\[ = —xj\fXjf — bXj\f_i. (3.38) 

The Toda model may be seen as an appropriate limit of the DST model (see also |53j). 
Indeed consider the following limiting process as 6 — 0: 

X„^e-«", x„^e''"(ri+p„) (3.39) 

It is clear that the harmonic oscillator algebra defined by (a:„, Xn, XnX^) reduces to the 
Euclidian Lie algebra (e^'^", Pn), and consequently the Lax operator reduces to: 



"'The associated Poissoii brackets for both DST and Toda models are defined as: 

dA dB dA dB 



( A ^ ( OA OB OA aB\ 



\dXn dXn dXn OXrJ ' 

'dA dB OA 

^dqn dpn dpn 



( , ] ^ f dA dB OA dB\ , , , , 
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where p„ are canonical variables. In this case the corresponding Hamiltonian may be 
readily extracted and takes the form 

TV 7V-1 
n=l n=l 

(2) 

and the corresponding An are expressed as 

In this case as well both formulas fl3.37p and (13.351) lead to the following set of equations of 
motions: 

Pn = Qn, Qn = e'^^+i""" - e''"-'?"-\ n G {2, . . . , iV - 1} 

Pn = Qn, Qn = e-'"'^ - q'^^-'^^-K (3.43) 



4 The continuous case 

4.1 Periodic boundary conditions 

Let us now recall the basic notions regarding the Lax pair and the zero curvature condition 
for a continuous integrable model following essentially [Hj. Define as being a solution of 
the following set of equations (see e.g. [6]) 

— = V{x,t,X)^ (4.1) 

— =Y{x,t,\)^ (4.2) 

U, V being in general n x n matrices with entries defined as functions of complex valued 
dynamical fields, their derivatives, and the spectral parameter A. Compatibility conditions 
of the two differential equation (14.11) . (14.21) lead to the zero curvature condition [3]-[5] 



U- V' 



U, V 



0, (4.3) 



giving rise to the corresponding classical equations of motion of the system under consider- 
ation. The monodromy matrix may be written from (14.11) as: 

T(x,?/,A) = Pexpj y l]{x',t,X)dx'y (4.4) 
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where apparently T{x,x,X) = 1. The fact that the monodromy matrix satisfies equation 
(14. ip is extensively used to get the relevant integrals of motion and the associated Lax pairs. 

Hamiltonian formulation of the equations of motion is available again under the r-matrix 
approach. In this picture the underlying classical algebra is manifestly analogous to the 
quantum case. Let us first recall this method for a general classical integrable system on 
the full line. The existence of the Poisson structure for U realized by the classical r-matrix, 
satisfying the classical Yang-Baxter equation (13. 5p . guarantees the integrability of the clas- 
sical system. Indeed assuming that the operator U satisfies the following ultralocal form of 
Poisson brackets 



T{x, y, A) satisfies: 

Ta{x,y,t,Xi), Tb{x,y,t,\2] 



rab(A-yu), Ua(x, A) +Ufe(2/,^) 6{x-y), 



rabi^i - A2), Ta{x,y,t,\i)Tb{x,y,t,X2) 

Making use of the latter equation one may readily show for a system on the full line: 

I In tr{T(x, 2/, Ai)}, \ntr{T{x,y, X2)}] = 

i.e. the system is integrable, and the charges in involution -local integrals of motion- are 
obtained by expansion of the generating function lntr{T[x, y, A)}, based essentially on the 
fact that T satisfies fl4.1D. 



(4.5) 



(4.6) 



(4.7) 



Let us now recall how one constructs the V-operator associated to given local integrals 
of motion. One easily proves the following identity using ( 14. 5p 

dM{x,X,fL) 



Ta{L,-L,X), lJb{x,fi) 
where we define 



dx 



M{x,L,-L,X,fi), Ub(x,/i) 



(4.8) 



M{x,X,^) =Tamx,X) rab{X-fi) Ta{x,L,X). (4.9) 

For more details on the proof of the formula above we refer the interested reader to [6] ; (14.80 
may be seen as the continuum version of relation (14. 6p . Recall that t{X) = trT{X) then it 
naturally follows from (14. 8p . and (14. 3p . that 

dY{x,X,fi) 



In t{X), U(x,A) 



dx 



+ 



V(x,A,/i), U(x,A) 



with 



Y{x,X,n) = t ^(A) traiTa{L,x,X) rab{X, fx) Ta{x,-L,X) 
and in the rational case (14. lip reduces to 

V(x, A, fi) = T{x, -L, A) r(L, x, A). 

A — ^ 



(4.10) 



(4.11) 



(4.12) 
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4.2 General integrable boundary conditions 

Our aim here is to consider integrable models on the interval with consistent "boundary 
conditions", and derive rigorously the Lax pairs associated to the entailed boundary local 
integrals of motion as a continuous extension of theorems 3.1.-3.3. For this purpose we 
follow the line of action described in [6], using now Sklyanin's formulation for the system 
on the interval or on the half line. We briefly describe this process below for any classical 
integrable system on the interval. In this case one constructs a modified 'monodromy' matrix 
T, based on Sklyanin's formulation and satisfying again the Poisson bracket algebras M or 
T. To construct the generating function of the integrals of motion one also needs c-number 
representations of the algebra M or T satisfying (13.14^ for SP and SNP respectively, such 
that: 

{k^{Xi), K±(A2)}=0. (4.13) 



Theorem 4-1-- The modified 'monodromy' matrices, realizing the corresponding algebras 



^ T, are given by the following expressions [11] (T being defined in fl3.16p ): 

T(x, y, t, A) = T{x, y, t, A) K'{X) f{x, y, t, A). 



(4.14) 



The generating function of the involutive quantities is defined as 

t{x,y,t,\)=tr{K+{X) r{x,y,t,X)}. 



Indeed one shows: 
Theorem 4-S.: 



|t(a;,?/,t, Ai), t(a;,?/,t, A2)| = 0, Ai, A2 G C. 



(4.15) 



(4.16) 



In the case of open boundary conditions, exactly as in the discrete integrable models, 
and taking into account (13.241) we prove 

(4.17) 



T,(0,-L,A), l]b{x,fi)j =M',{x,X,^)+ ^Ma{x,X,fi), Ub(x, ^) 
where we define 

M(x,A,/i) = T{0,x,X)rab{X- fi)T{x,-L,X)K^{X)f{0,-L,X) 
+ no, -L, X)K~{X)f{x, -L, A)f:,(A + ^)r(0, x, A). 



(4.1^ 
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Finally bearing in mind the definition of t{X) and fl4.17p we conclude with: 



Theorem 4-3.: 



where 



dY(x X u) 

In t(A), U(x, fi) \ = + A, /i), U(x, //) (4.19) 



Yix,X,fi)=t-\X) trjK+iX) Maix,X,fi)). (4.20) 



As in the discrete case particular attention should be paid to the boundary points x = 
0, — L. Indeed, for these two points one has to simply take into account that T{x, x, X) = 
T{x,x,X) = I. Moreover, the expressions derived in fl4.18p . fl4.20p are universal, that is 
independent of the choice of model. As was remarked upon when discussing the discrete 
case, a quadratic algebra of the form f l3.14p was initially obtained in the continuous case 
[15] when extending the derivation of 4.1 to situations where the Poisson brackets (14. 5 p are 
non-ultralocal, exhibiting 6'{x — y) terms. Connection to "boundary" effects was discussed 
previously (see Section 3.2). 



4.3 Example 

We shall now examine a particular example associated to the rational r- matrix (12. 7p . that 
is the glj\j- NLS model. Although in [13] an extensive analysis for both types of boundary 
conditions is presented, here we shall focus on the simplest diagonal {K^ = I) boundary 
conditions. The Lax pair is given by the following expressions [6], [51] : 



U = Uo + AUi, V = Vo + AVi + A^Vs (4.21 



where 



1=1 i=l 

Vo = « ^ (ipiipjEij - lipil"^ E^f^) - {^{ip'^EiX - ipiEf^i), 

i, j=l i=l 

Vi = -Uo, V2 = -Ui (4.22) 
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and ipi, ipj satisfja: 

^j{y)^ = [^p^{x), = 0, {Mx), ^Pj{y)^ = 6,.j 5{x-y). (4.24) 

Note that we have suppressed the constant k from fl4.22p compared e.g. to [I3] by rescahng 
the fields (^/^i, v^(^i, ipi). 

From the zero curvature condition (14. 3 p the classical equations of motion for the gener- 
alized NLS model with periodic boundary conditions are entailed i.e. 

^djHi^ = -^^^If^ + 2KY,m^,t)\'U^,t), J G {1, . . . ,Ar- 1}. (4.25) 

3 

It is clear that for A/" = 2 the equations of motion of the usual NLS model are recovered. 
The boundary Hamiltonian for the generalized NLS model may be expressed as 

n= 5^(Kiv^.(x)p5^i^,(x)r+^:(x)^:(x) 

i=i j=i 

- E (v':(o)^.(o)+^.(o)^:(o)) + J2 (i^[{-m{-L)+M-m{-L)).{4.2Q) 



i=l i=l 



One sees here that the i^"- matrix indeed contributes as a genuine boundary effect. The 
Hamiltonian, obtained as one of the charges in involution (see e.g. [13] for further details) 
provides the classical equations of motion by virtue of: 

^^^ = {7^(0,-L), Ux,t)}, ^M^ = |7^(o,-L), ^,{x,t)], 
-L<x<0. (4.27) 

Indeed considering the Hamiltonian Ti, we end up with the following set of equations with 
Dirichlet type boundary conditions 



.dilJi{x,t) d^i)i{x,t) 



2 , / ,N A^-l 



l- 



dt d^x 



2/cElV',-(x,t)pV'i(2;,t) 



^,(0) = 7/;i(-L) = {!,..., AT- 1}. (4.28) 

For a detailed and quite exhaustive analysis of the various integrable boundary conditions 
of the NLS model see [13]. Note also that the M = 1 case was investigated classically on the 
half line in [55], whereas the NLS equation on the interval was studied in |56j . 
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The Poisson structure for the generahzed NLS model is defined as: 

L 



r _ -s:^ i ^ ( 6 A SB ^ 
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As mentioned our ultimate goal here is to derive the boundary Lax pair, in particular the 
V operator. In general for any glj^ r- matrix we may express fl4.20p . taking also into account 
f l4.19p . for the two types of boundary conditions already described in the first section, i.e. 

W(x,X,^) = —^T(x,-L,X)K~(X)f(0,-L,X)K+(X)T(0,x,X) 
A — /i 

+ ^-^r(0,x,A)K+(A) r(0,-L,A)/\-(A)r(x,-L,A) forSP (4.29) 
A + /i 

and for the SNP boundary conditions we obtain: 

V(x,A,^) = — ^T(x,-L, A)ir-(A)f (0,-L, A)K+(A)T(0,x, A) 
A — /i 

+ ^^^T*(x, -L, X)K~\X) T*(0, -L, X)K+\X)f\0, x, A) for SNI?4.30) 
A + /i 

Again for the boundary points x = 0, —L we should bear in mind that T{x, x, A) = 
T{x,x,X) = I. Ultimately we wish to expand T(A), T(A) in powers of A~^ in order to 
determine the Lax pair for each one of the integrals of motion. For a detailed description of 
the derivation of the boundary integrals of motion for the generalized NLS models see [13] . 
Hereafter we shall focus on the SP case with the simplest boundary conditions i.e. = I. 
Expanding the expression fl4.29l) in powers of A~^ (we refer the interested reader to Appendix 
C for technical details) we conclude that Y^^\x, A) -the bulk part- coincides with V defined 
in fl4.2ip . fl4.22p . and for the boundary points xi, G {0, —L} in particular: 

2 Af-l JV-1 Af-1 



=1 i,j=l i,j 



V^^^xt, A) = -— ( ^ E,, - Ej^^^ + i iPi{xh)iJjj{xb)Eij - % ^ {^\{xb)EiM - iJJi{xh)Ex,^ . 



(4.31) 



We may alternatively rewrite the latter formula as: 

Y^^\xb, A) = Y{xb, X)+tJ2 \i^i{xb)\^Exx + Xj2iMxb)Ei^ + iJi{xb)EMi). (4.32) 

1=1 1=1 

The last two terms additional to V (14. 2 ID . fl4.22p are due to the non-trivial boundary condi- 
tions; of course more complicated boundary conditions would lead to more intricate modifi- 
cations of the Lax pair V, however such an exhaustive analysis is beyond the intended scope 
of the present investigation. It can be shown that the modified Lax pair (U, Y^'^^) gives 
rise to the classical equations of motion fl4.28p . It is clear that the 'bulk' quantity V in the 
case of SP boundary conditions remains intact. In the SNP case on the other hand we may 
see that even the bulk part of the Lax pair is drastically modified, due to the fact that the 
bulk part of the corresponding integrals of motions is also dramatically altered. We shall 
not further comment on this point, which will be anyway treated in full detail elsewhere. 
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A Appendix 

We present here technical details on the derivation of the conserved quantities for the gener- 
alized NLS model on the full line (see also [13]). Recall that for A — > ±zoo one may express 
T as [6] 

T{x, y, A) = (I + W{x, A)) exp[Z{x, y, A)] (I + W{y, X))-' (A.l) 

where W is an off diagonal matrix i.e. W = Yli^j ^ijEij, and Z is purely diagonal Z = 
Y^-^iZiiEii. Also 

oo „(n) oo TJ/W 

n=—l n=l 

The first step is to insert the ansatz flA.ip in equation (14.11) . Then we separate the diagonal 
and off diagonal part and obtain the following expressions: 

Z' = AUi + {VoWY^'^ 

W' + WZ' = l]o + (VoWY^^ + XViW (A.3) 

where the superscripts (D), (O) denote the diagonal and off diagonal part of the product 
UoVF. Recall that W = ^ij^ij, Z = J2i ^uEu then it is straightforward to obtain: 

i=l 

(UoW^)(^) = {^Wu^Ei, + ^iW^.Eu^ . (A.4) 

Substituting the latter expressions (1A.4I1 in (1A.3I) . we obtain 

Z{L, -L, A) = -zAL(^ - ^^^) + / dx(iP,W^iEu + ipiWi^E^^^^ . (A.5) 

i=i i=i 

The leading contribution in the expansion of (In trT), (where T is given in (lA.ip ) for i\ ^ oo 
comes from Zj^j^, with a leading term iXL (all other Zu, i ^ M have a —iXL leading term, 
so when exponentiating such contributions vanish as iX oo), indeed 

Zxx{L,-L,X) = iXL + y/^Y dx iJji{x)Wix{x). (A.6) 
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Due to (1A.6P it is obvious that in this case it is sufficient to derive the coefficients Wij\f. In 
any case one can show that the coefficients Wij satisfy the following equations: 



J2 -^^Yl (^AT.^AT. - W,^fE,^') +V^Y1 {^^W^^E^f^ + i^^WfJ^E, 

[i^jWMjWijEij + ^|JiWi^Wj^fEj^f) ■ 



(A.7) 



(n) 



Finally setting Wij = Xl^i ^'^d using (1A.7I) we find expressions for Wl]:^ i.e. 



X) 



In the boundary case we shall need in addition to flA.SP the following objects: 



W 



(1) 



(2) 



-iW 



'(1) 



i^h ij^N, jjtN 



(A.9) 



Based on the latter formulas, the expression ( lA.ll) . and defining A) = W{x,—\) we 

may rewrite fl4.29p as: 

^ '1 + + w{x)Y^ + -^(1 + w{x))Ej^m{^ + 



A — /i A + /i 

V(0,A,/.) = W)-^(-^X+K+(A) + -^K+(A)X 

V(0, A,^) = {X^)-^{^K-{\)X- + ^X-K-(A) 



A — /i 



X — jj 



(A.IO) 



where we define: 



X+ = (1 + t^(0))E^^(l + iy(0))-\ = (1 + W^(-L))SAr^(l + W{-L)y' 



[1 + W{0, A))-^K+(A)(1 + H^(0, A)) 



MM 



[1 + Vr(-L, A))-^A-(A)(1 + W{-L, A)) 



(A.ll) 
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